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ABSTRACT 

We present combined Submillimeter- Array (SMA) + single-dish images of the (sub)millimctcr dust 
continuum emission toward two prestellar cores SMI and B2-N5 in the nearest star cluster forming 
region, p Ophiuchus. Our combined images indicate that SMI and B2-N5 consist of three and four 
condensations, respectively, with masses of 1CP 2 — 10 -1 M© and sizes of a few hundred AU. The 
individual condensations have mean densities of 10 8 — 10 9 cm~ 3 and the masses are comparable to 
or larger than the critical Bonner-Ebert mass, indicating that the self-gravity plays an important role 
in the dynamical evolution of the condensations. The coalescence timescale of these condensations 
is estimated to be about 10 4 yr, which is comparable to the local gravitational collapse timescale, 
suggesting that merging of the condensations, instead of accretion, plays an essential role in the star 
formation process. These results challenge the standard theory of star formation, where a single, 
rather featureless prestellar core collapses to form at most a couple of condensations, each of which 
potentially evolves into a protostar that is surrounded by a rotating disk where planets are created. 
Subject headings: ISM: clouds — ISM: kinematics and dynamics — ISM: structure — stars: formation 



1. INTRODUCTION 

According to the standard theory of star formation, 
stars form o ut of dense cor e s embedded in molecula r 
clouds (e.g.. IShu et all [19871: IMcKee fc Os"trikerl [2007h . 
The dense cores prior to the protostellar formation, or 
prestellar cores, ar e thought to be smoo th with no con- 
densations inside ([Bacmann et al.1 120 00). The gravita- 
tional collapse of such a core impedes the growth of the 
local density fluctuations, resulting in the formation of 
at most a couple of condensations, each of which poten- 
tially evolves into a protostar that is surrounded by a 
rotating disk where planets are created. Recent millime- 
ter and submillimeter observations have revealed that the 
core mass functions (C MFs) resemble the stellar initial 
mass function (IMF) dMotte et all 119981: iStanke et all 
120061: iMaruta et al.ll2010D . The resemblance between the 
CMFs and IMF appears to claim that these cores are 
the direct precursors of individual stars or small stellar 
systems like binaries. 

The above scenario of star formation implicitly as- 
sumes that the core has a smooth density distribution 
and does not contain significant condensations. This 
assumpt ion appears to be support ed by recent observa- 
tions by iSchnee et al.l (|2010l [20121 ). who suggested that 
the prestellar cores do not contain significant substruc- 
tures on the basis of the interfero metric observation s 
with arcsecond re solution (see also lOlmi et ail l2005f) . 
iMaurv et al.l (|2010T ) studied substructures around several 
Class sources with subarcsecond resolution and found 
no clear signs of fragmentation at 100 AU. On the other 
hand, several previous observations revealed the pres- 
ence of gravitationally-unbound, substellar-mass struc- 
tures with sizes of a few thousands AU inside prestel- 
lar cores in Taurus (lLanger et al.lll995l : lPeng et al.lll998l : 
iTakakuwa et al.ll2003HRov et aLlboTl . The origins and 
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fates of these small structures inside the prestellar cores 
remain to be elucidated. These previous studies of sub- 
structures in prestellar cores have concentrated primarily 
on the relatively-isolated environments. Hence, none of 
the previous observations have clearly answered the ques- 
tion as to whether prestellar cores in cluster-forming re- 
gions have significant substructure or not, because of the 
lack of high spat i al reso lution observations [see, however, 
iKamazaki et al.l (|2001[ ) for the discovery of thousand- 
AU scale, substellar-mass structures in Oph A]. However, 
since most stars form in the clustered environments, it is 
important to elucidate the structures and properties of 
prestellar cores in the clustered environments. 

In the clustered environments, the identified cores tend 
to be more compact than those in quiescent regions, and 
their sizes are comparable to a telescope beam size. An 
area covered by a typical single prestellar core in quies- 
cent regions is often occup ied by a small clus ter of cores 
in cluster-forming regions (|Motte et al.l 11998) . To verify 
whether prestellar cores in the cluster environments con- 
tain substantial internal structures or not, we analyze in 
this Letter two representative prestellar cores (SMI and 
B2-N5) in p Ophiuchus, the nearest cluster-forming re- 
gion at a distance of 120 pc, using the dust continuum 
data obtained by SMA. The SMA data have spatial res- 
olution of a few hundreds AU, one of the finest among 
the available data of prestellar cores. 

The first target SMI, located in the middle of the Oph 
A filamentary cloud, has the strongest submillim eter con- 
tinuum emission in p Oph (jAndre et al.lll993l) . The fil- 
amentary morphology is suggestive of interaction with a 
nearby young B star SI. The total mass and dust temper- 
ature of SMI are estimated to be 2 Mq and Td m 20 K, 
respectively (|Motte et al.l 119981 ). The second target B2- 
N5 is located in the Oph B2 cloud that contains a number 
of prestellar cores. In the Oph B 2 cloud, a gigantic out - 
flow has recently been detected (Nakamu ra et al.ll2011l) . 
This core was identified in the N2H + (1-0) emission and 
is one of the cores having the largest column densities 
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in the Oph B2 cloud. The total mass an d temperature 
are e stimated to be about a fe w tenth M<^ ( Friesen et al.l 
[20101 N 2 H+) and T « 15 K (|Friesen et all 120091 NHA 
respectively. 

The rest of the paper is organized as follows. The de- 
tails of our data are described in Section [2] The results 
and discussion are presented in Sections [3] and 01 respec- 
tively. 

2. DATA 
2.1. Single-dish Continuum Data 

The p Oph cluster-forming region consists of several 
sub-parsec scale clouds with masses of order of 10 2 Af Q , 
each of which contains several prestellar cores. Our tar- 
get cores SMI and B2-N5 are located in the Oph A and 
B2 clouds, respectively. Here, wc briefly describe the 
dust continuum data toward p Oph, taken by the single- 
dish observations. 

The 850 pm data of the Oph A cloud taken with the 
SCUBA were obtained through the archiv al system of the 
COMPLETE Survey (IRidee et al.l [2006h . The FWHM 
beam size of the JCMT was 14" at 850/im. The mean 
rms noise level of the data is ~10 mJy beam -1 . 

The 1.1 mm data of the Oph B2 cloud were obtained 
with the AzTE C camera dWilson et al.l 120081) on the 
ASTE telescope (jEzawa et al 112004ft . The FWHM beam 
size of the ASTE was 28" at 1.1 mm. The observa- 
tions were performed in the raster scan mode. Each 
field was observed several times with azimuth and ele- 
vation scans. The separation among scans was adopted 
to be 117", which is a quarter of the AzTEC field of 
view (FoV ~ 7'. 8). The effective beam size of the data 
was 40" after the FRUIT imaging which is an iterative 
mappin g method to reco ver the spatially-extended com- 
ponent ( Liu et all [2010). The noise level was less than 
10 mJy beam -1 in the entire observed area. The details 
of the data will be presented elsewhere. 

In general, it is difficult to completely recover the emis- 
sion from extended structures with bolometric observa- 
tions using the ground-based telescopes due to the at- 
mospheric emission. For SCUBA, the emission from the 
structures larger than ~ 2' is suppressed in removing the 
atmospheric emission during the data reduction process, 
making the 850 pm map mostly devoid of the extended 
emission. The AzTEC image also has the same problem 
for the structures larger than the field of view of ~ 8'. 
However, this effect is expected not to be severe because 
the sizes of the cores and the condensations identified be- 
low are smaller than - 2' and 8' for SCUBA and AzTEC, 
respectively. 

2.2. SMA Data 

The target prestellar cores SMI and B2-N5 were ob- 
served in the SMA compact configuration at 870 pm and 
1.1 mm, respectively. We obtained the continuum data 
through the SMA da ta archive. The SMA is a double- 
sideband instrument ()Ho et al.l I2004D . having a number 
of spectral windows of the SMA correlator ("chunks"). 
Each target was observed in a single pointing mode. 

Table [1] summarizes the observational parameters. The 
minimum projected baseline lengths were 11 and 12 fcA 
(8.2" and 7.5") at 870 pm and 1.1 mm, respectively, and 
for a Gaussian emission distribution with a FWHM of 



~ 7" (~ 800 AU), the peak flu x recovered is ~ 50% o f 
the peak flux of the Gaussian (|Wilner fc Welch! 11994ft . 
This problem is called the missing flux problem. To re- 
cover the missing fluxes from extended structure and ob- 
tain more reliable images and physical properties of the 
structures, we combin ed the SMA data and the avail- 
able single-dish data (Takakuwa ct al. 2007): i.e., 850 
pm SCUBA data for SMI and 1.1 mm AzTEC data 
for B2-N5. Below, we mainly discuss the combined 
SMA+single-dish images. 

3. RESULTS 

The combined images are presented in Figure [T] For 
comparison, the original SMA images are also shown in 
Fig. [2j Our high spatial dynamic range images con- 
structed by combining the single-dish data and SMA 
data have spatial resolution and sensitivity to dust emis- 
sion superior to those of existing millimeter and sub- 
millimeter observations toward prestellar cores in active 
cluster-forming regions, without significant effects of the 
missing flux inherent in any interferometric observations. 
Both images have clearly revealed the presence of signif- 
icant substructure inside the prestellar cores. Here, we 
define a compact 10 2 — 10 3 AU-scale structure with a peak 
intensity larger than 10 a rms noise level as a condensa- 
tion. For SMI, we identified three condensations, labeled 
as al, a2, and a3 in Fig. QJa). The southernmost conden- 
sation, al, appears extremely compact, and the internal 
structure is unresolved even with a few arcsecond resolu- 
tion. The existence of the thousands AU-scale substruc- 
ture in the Oph A cloud has been poi nted out by the 
previ ous interferometric observations fKamazak i et all 
12001ft . Our image indicates that the substructure pre- 
viously identified as the stro ngest 3mm continuum e mis- 
sion (core A in Figure 2 of Kamaz aki et al.l (]2001ft ) is 
much more compact and consists of two condensations 
al and a2. The substructure is more prominent in B2- 
N5, which contains several small condensations, as shown 
in Fig. [lib). Here, we identified 4 condensations in the 
map. The separations of the condensations are about 
5—10 arcseconds, or 10 3 AU. 

We note that the condensations identified from the 
combined images can be recognized easily even in the in- 
terferometric images alone (Figure^, and the small-scale 
structures are essentially the same as those presented in 
Figure [T] The peak fluxes of all the condensations iden- 
tified in the SMA images alone are more than 4<r rms 
noise levels, indicating that these condensations are the 
real structures existing in the prestellar cores. 

Table [2] summarizes the physical properties of the iden- 
tified condensations. We applied two-dimensional Gaus- 
sian fitting to the condensations by using the IMFIT task 
in MIRIAD and derived their sizes, positions, and peak 
and total fluxes. The masses of the condensations are 
estimated using the dust opacity of K\ = . 1(250 /xm/A)^ 
g -1 cm 2 with = 2 (jHildebrand et al.l Il98l . This 
opacity is in agreement with the values adopted by 
other stu dies of prestellar cores within a factor of a 
few (e.g..lOssenkopf fc Hennind[l99l lEvans et al.ll2009l: 
lAndre et al.lll996ft . The typical masses of the condensa- 
tions are of the order of 10 -2 — 10 -1 M Q . For compari- 
son, the mass of the critical Bonner-Ebert sphere whose 
mean density coincides with the estimated mean density 
of each condensation is listed in Table [2] The elongated 
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beam shapes influence the estimation of the local den- 
sities and therefore Mbe- Since the identified conden- 
sations have sizes comparable to the beam sizes, Mbe 
presented in Table [2] may be underestimated. Almost all 
the condensations have masses comparable to or larger 
than the critical Bonner-Ebert masses (Mbe), indicat- 
ing that the self-gravity plays an important role in the 
dynamical evolution of the condensations. 

We note that the masses estimated from the continuum 
emission depend strongly on the adopted opacity and 
temperature. These two parameters are generally diffi- 
cult to measure accurately. The dust opacity has uncer- 
tainty at least by a factor of a few. Therefore, the masses 
of the condensations listed in Table [2] are reduced by a 
factor of a few when the actual dust opacity is smaller. 
The adopted temperatures also have uncertainty. Here, 
T = 20 K is adopted for SMI , although a large r tem- 
perature of 27 K is reported (|Andr£eFlil| [T993) . For 
B2-N5, the dust temperature is assumed to be equal to 
the kinetic temperature o f 15 K, which is d etermined 
by the NH 3 observations ([Friesen et al.l I2009T). In Oph 
B2, t he protostellar outflow from EL32 (jNakamura et all 
1201 If ) may have directly injected turbulent motions in 
adjacent dense regions including B2-N5, and may have 
increased the temperature of B2-N5. However, the N2H + 
(J = 1 — 0) hyperfine fitting indicates that the t emper- 
ature of B2-N5 appears low (|Friesen et all 120091 ) . The 
NH3 observations may not trace the densest gas in Oph 
B2, and the adopted temperature may be overestimated 
in the densest parts. Even if we take into account these 
uncertainties of the dust opacity and temperature, all the 
condensations still have masses that are comparable to or 
larger than Mbe, for SMI. For Oph B2, the masses of the 
condensations may become somewhat smaller than Mbe 
when we take into account the larger dust opacity and 
smaller temperature. Even in that case, the existence of 
the large-amplitude fluctuations at this small-scale is sur- 
prising because such large-amplitude fluctuations give a 
significant impact on the structure form ation during the 
gravi tational collapse of the parent core ([Goodwin et all 
12001 . 

4. DISCUSSION 

The significant substructure inside the prestellar cores 
implies that the basic unit of star formation in the clus- 
tered environments is much smaller than previously con- 
sidered. In fact, recent theoretical and observational 
studies suggest that prestellar cores formed in the clus- 
tered environments tend not to be self-gravitating, but 
pressu re-confined because of the strong inter-cor e turbu- 
lence (|Maruta et alll2O10HNakamura fc Lill201lD . 

How were these substellar-mass condensations created 
in parent prestellar cores? One possibility is the fragmen- 
tation due to internal turbulent motions, or turbulent 
fragmentation, in short. In the standard theory of star 
formation, the cloud turbulence is considered to be less 
important inside the prestellar cores. However, cluster- 
forming regions are highly turbulent. Supersonic tur- 
bulence potentially triggers the formation of dense re- 
gions where the turbulent d issipation preferentially hap- 
pens (jNakamu ra fc Lill201lD . In such dense regions the 
large- amplitude fluctuations can be left. To quantify the 
dynamical importance of the core internal turbulence, 
we present in Figure [3] the ratios of the turbulent to 



gravitational energy (aturb = U/W) for the prestellar 
cores in p Oph as a function of the LTE mass, MltEj 
where U is the turbulent energy U = MuteSv^ st /2, W 
is the gravitational energy W = GM 2 TE /i? corc , Sv^t 
is the three-dimensional nonthermal velocity dispersion, 
G is the gravitational constant, and i? C orc is the core 
radius. All the prestellar cores were identified in a 
high density tracer, H 13 CO+ (J = 1-0), using the 
Nob eyama 45m telescope whose beam size is about 2000 
AU ([Maruta et al.1 l2010h . The mean value of aturb is 
estimated to be about 0.8, suggesting that the internal 
turbulence can contribute significantly to the core dy- 
namics. This is in contrast to the core s of quiescent 
regions like Taurus ([Goodwin et al.l 120041 ) . where a tU rb 
tends to be smaller. In addition, the values of aturb are 
significantly larger than th e threshold for turbul ent frag- 
mentation Of aturb >0.03 (|Goodwin et alj|200l . Thus, 
the turbulent motions are responsible for the formation 
of the condensations inside the prestellar cores in p Oph. 

H 13 CO + sometimes suffers from the molecular de- 
pletio n in the dense pa rts of the prestellar cores. 
Since iMaruta et all (|2010D adopted a constant H 13 CO+ 
fractional abundance, the estimated core masses (and 
thus a tU rb) may be underestimated (overestimated) if 
H 13 CO + is significantly depleted. However, the deple- 
ti on effect is exp ected to be minor as discussed in §4.3 
of IMaruta etafl (|2010h . Therefore, we believe that the 
estimation of the core masses and aturb is reasonable, 
and our conclusion is not altered by the effects of the 
molecular depletion. 

What is the fate of these substellar-mass condensa- 
tions? Both the prestellar cores contain condensations 
with sizes of order 10 2 AU and masses of a few hun- 
dredth M Q . Almost all the condensations have masses 
comparable to or larger than the critical Bonner-Ebert 
mass. Therefore, they are expected to be gravitationally 
bound. From the mean separation between the conden- 
sations of a few xlO 2 AU the coalescence timescale of 
these condensations is estimated to about 10 yr, compa- 
rable to the free-fall times of the individual condensations 
with densities of 10 s cm -3 , where the relative velocities 
among condensations are assumed to be equal to the ve- 
locity dispersion of B2-N5, Sv ~ 0.2 km s _1 . Therefore, 
these condensations can either merge with themselves or 
collapse individually to form small groups of low-mass 
protostars and/or prot-brown dwarfs. Here, the mean 
separation projected on the plane of the sky is adopted 
to calculate the coalescence timescale and therefore the 
estimated coalescence timescale gives a lower limit. Un- 
less the actual mean separation is larger than the adopted 
value by a factor of ~ 10, the coalescence timescale is still 
comparable to the local free-fall time and thus our con - 
clus ion is not altered signifi cantly (see lPeng et all (jl998f ) 
and lTakakuwa et all ([2003) for more quantitative analy- 
sis on the separations of condensations.). 

According to the standard theory of star formation, a 
typical prestellar core has sizes of 10 4 AU and masses 
of a few M Q — 10 M and is considered to have no 
significant substructures. The fragmentation during the 
prestellar phase is considered to be significantly impeded 
by the glob al gravitational collapse before a opaque disk 
is for med (|Hennebelle fc Tevssierl 120081 : Machid a et all 
2008]). The prestellar cores could fragment into at most 
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a couple of condensations that grow in mass by accre- 
tion and eventually evolve into protostars. Thus, a bi- 
nary or small multiple system is created. The rotat- 
ing disks are created around individual protostars and 
the planets are produced there. However, the recent 
hydrodynamic simulation of gravitational collapse of a 
prestellar core has posed a challenge to this standard 
picture. In a series of three-dimensional radiation hydro- 
dynamic simulations, the protostellar radiation destroys 
the o paque disks that would evolv e into circumstellar 
disks (jSchonke feTscharnuterll2011[ ). Although whether 
the disks are completely destroyed or not may depends 
on the size of the disk formed, this effect of the proto- 
stellar radiation significantly constrains the planet for- 
mation from the circumstellar disks. In our scenario, 
substellar-mass condensations are already formed before 
the protostars are created. Some of these substellar-mass 
condensations potentially evolve into giant planets orbit- 
ing around protostars or proto-brown dwarfs under the 
influence of the protostellar radiation that can remove a 
significant amount of mass from the condensations. 

Here, we revealed that the two representative prestel- 
lar cores in the p Oph cluster forming region contain 
hundreds AU-scale structures. Another compact struc- 
ture with ^1000 AU was also found toward the southern 
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TABLE 1 

Parameters for the SMA Observations of SMI and B2-N5 



Parameter 


SMI 


B2-N5 


„, „ 

W DSSrvlIlg DcitrG 


July zy, zuu ( 


Tnno 1 900& 
j Line r , zuuo 


NuiiibGr of ctiitrcniicis 


7 


e 
o 


R A ( 12000^1 


16-26-27 60 


16-27-29 15 


Decl.(J2000) 


-24:23:55.0 


-24:27:02.0 


Primary beam HPBW 


36" 


45" 


Synthesized beam HPBW 


2.7" x 1.3" (P.A.=49°) 


2.7" x 2.2" (P.A.=40°) 


Baseline coverage 


10.9 - 148.2 (kA) 


12.1 - 114.2 (kA) 


LO frequency 


341.5 GHz 


274.2 GHz 


Bandwidth 


4.0 GHz 


3.8 GHz 


Gain calibrator 


NRAO530 


1625-254 


Flux of the gain calibrator 


1.4 Jy 


0.79 Jy 


Passband calibrator 


3C273, 3C454.3 


3C279, Jupiter 


Flux calibrator 


Uranus 


Titan 


System temperature 


150 - 600 (K) 


300 - 600 (K) 


rms noise level 


4.8 mjy/beam 


1.2 mjy/beam 



TABLE 2 

Physical properties of condensations in prestellar cores. 



Name 


R.A. 


Dec. 




St 


Size 


Mass 


density 






(J2000) 


(J2000) 


(mjy/beam) 


(mjy) 


(AUxAU) 


(Mq) 


(xlO 8 cm" 3 ) 


(M e ) 


al 


16:26:27.8 


-24:24:00 


455 


742 


350 x 167 


0.15 


30.7 


0.014 


a2 


16:26:27.7 


-24:23:57 


176 


958 


770 x 287 


0.17 


4.6 


0.036 


a3 


16:26:26.6 


-24:23:54 


138 


769 


591 x 532 


0.13 


2.2 


0.053 


bl 


16:27:29.6 


-24:27:07 


15.2 


44 


532 x 310 


0.033 


1.4 


0.03 


b2 


16:27:29.5 


-24:27:04 


17.4 


136 


960 x 702 


0.10 


0.53 


0.056 


b3 


16:27:28.9 


-24:27:02 


15.8 


47 


548 x 340 


0.035 


1.3 


0.036 


b4 


16:27:28.6 


-24:26:59 


13.1 


49 


541 x 480 


0.036 


0.79 


0.046 



Note. — The physical quantities listed in the table are evaluated from deconvolved Gaussian-fitted model parameters. The core masses are 
estimated by adopting the dust temperatures of 20 K and 15 K for the prestellar cores, SMI and B2-N5, respectively. The size is given by the 
FWHM major and minor axes obtained from the IMFIT Gaussian fitting. 



6 



i.rstr 



24' U 



-24 "Mitt* 




r 


"■ ■- — i ■ 


J 


| IT lJ 


i a.i 


ft 


I" 









2r«- 


(t) Oph 




M 












1 M-i- 


i 







H,iii*t ItnJIh Unite liUMQi 




ITra.lfi* J7m5fc IThH* IbhllaiiOi 



Fig. 1.— (a) Combined SCUBA and SMA image toward SMI . The contour levels are -20, 60, 80, 100, 130, 160, and 320 mjy/beam. 
The cross is the position of SMI identified by Mottc ct al. (1998). The alphabets designate the identified condensations, (b) Combined 
AzTEC/ASTE and S MA image toward B 2-N5. The contour levels are -3, 8, 10, 12, 14, and 16 mjy/beam. T he cross is the positio n 
of B2-N5 identified by Fricscn ct al. (2010). (c) the SCUBA 850/im continuum image toward the Oph A region ( Johnstone ct al. 120001 ). 
The white and black filled circles are the B star, SI, and the prototypical Clas s YSO, VLA 1623, respectively. The positions of some 
submillimeter continuum sources are indicated by the squares (Mottc ct al. 1998). (d) the AzTEC/ASTE 1.1 mm continuum image toward 
the Oph B2 region. The black filled circles are the positions of Class I YSOs, EL32 and EL33. In panels (a) and (b), the synthesized beams 
are shown in the lower right of the panels. 
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Fig. 2. — (a) Original SMA image toward SMI. The contour levels are -20, 40, 60, 80, 100, and 200 mjy/beam. The cross is the position 
of SMI. The alphabets designate the condensations identified in Fig. [TJa). (b) Original SMA image toward B2-N5. The contour levels are 
-3, 4, 6, 8, 10, 12 mjy/beam. The alphabets designate the condensations identified in Fig. [l^b). The synthesized beam is shown in the 
lower right of each panel. 
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Fig. 3. — The ratio of turbulent to gravitational energy («turb) as a function of LTE mass (A/lte) for the p Oph prestellar cores identified 
from H 13 CO+ (1-0) UMaruta et alj|201Ch . 



